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Abstract 

O-h In a very recent paper "On eigenvalues and equivalent transformation of trigonometric matrices" (D. 

Zhang, Z. Lin, and Y. Liu, LAA 436, 71-78 (2012)), the authors motivated and discussed a trigonometric 
matrix that arises in the design of finite impulse response (FIR) digital filters. The eigenvalues of this 
matrix shed light on the FIR filter design, so obtaining them in closed form was investigated. Zhang et 
al. proved that their matrix had rank-4 and they conjectured closed form expressions for its eigenvalues, 
leaving a rigorous proof as an open problem. This paper studies trigonometric matrices significantly more 
general than theirs, deduces their rank, and derives closed-forms for their eigenvalues. As a corollary, it 
yields a short proof of the conjectures in the aforementioned paper. 

1 Introduction 

Matrices generated by trigonometric functions enjoy a wide-range of appealing analytic and numerical prop- 
erties, whereby they are important ac ross a wide yariety of problems. The present paper is motivated by the 
very recent article (in this journal) of IZhang et al. 2012j . who studied certain trigonometric matrices arising 



in the design finite impulse response (FIR) digital filters. In that paper, the authors motivated the importance 
of obtaining ei genvalues of these m atrices for their FIR design task. 

Specifically, IZhang et al. 2012| considered trigonometric matrices of the form 



PH := 



A(u) B(w) 
B(co) T A{u) 



(1.1) 



where A(ui) and B{uS) are n x n matrices (n > 2), and w is the digital frequency variable with < U) < 2n. 
The entries of A{uj) and B(oj) (denoted with lowercase letters) are given by 

a%j{(jS) := h(i + j — 2) cos(io; — jui), and bij(uj) := h(i + j — 2) sin(io; — joS). (1-2) 

For the matrix P(u>) defined using (fl~2l) . IZhang et all (20121 ] made the following conjecture. 



Conjecture 1 ( [Zhang et aTll2012j ). P(u) has one positive eigenvalue A+ and one negative eigenvalue A_, 
both of which are of multiplicity 2, are independent of uj, and are given by 



K = - A [n-l + ] / § j A_ = i (n-l-y j. (1.3) 

The other eigenvalues of P(u>) are zeros. 

In this paper, we consider trigonometric matrices more general than (II. 2j) . and derive their eigenvalues in 
closed form. As a corollary (Corollary [S]) we obtain a shor t proof of Conjectu re ffl Our more general matrices 
may be useful in applications such as those mentioned bv lZhang et al. 2012| and references therein. 
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2 Generalized trigonometric matrices 

In this paper we consider generalized trigonometric matrices 



P(u) 



A(u) B(u) 
B t (uj) A(u) 



(2.1) 



that have the same form as (11.11) but with more general component matrices A(u>) and B(u>) defined by 

Oij(w) := hj cos(xi — Xj), and bijiw) := Uj sin(xi — Xj), (2-2) 



where Xi and Xj are components of an arbitrary vector x € R™, and ijj are the entries of an arbitrary, symmetric 
rank-2 matrix L defined as 

L:=lh T + hl T , for vectors I, h e M™. (2.3) 

In particular, matrix (jl.ip corresponds to L = Ze T + e/ T for vector I having entries ij = (i — l)/2, where e 
denotes the all ones vector. 

To analyze the rank of P(uS) and its eigenvalues we need to introduce more notation. First, we drop ui, 
because the end results do not depend on it; then, we introduce the pure trigonometric matrix 



P := 



A B 
B T A 



where the entries of the component matrices A and B are given by 

ciij := cos(xi — Xj), and bij := sin^ — Xj), 
Using this notation we see that the matrix (|2.1I) can be written as 



P = PQ 



1 1 
1 1 



L ) =P(D 



L L 
L L 



(2.4) 



(2.5) 



(2.6) 



where €5 is the Kronecker and the Hadamard product. But we know that rank(X 7) < rank(X) rank(F) 
for any X and Y; moreover rank([^ ^1) = rank(i) = 2, whereby from (|2.6p it follows that 



rank(P) < 2rank(P). 



(2.7) 



We note that inequality (|2.7j) may be tightened to the equality rank(P) = rank(L) rank(P) (a fact that follows 
easily once we have derived the eigenvalues of P). But for now, the inequality suffices and we proceed to 
analyze the rank of P. To that end, Lemma [5] proves useful. 



Lemma 2. Let A = UU T , where U is any n x 2r (2r < n) matrix of rank 2r. Let J^r be the 'symplectic 
identity' 

' I r ' 
-Ir 



where I r is the r x r identity matrix. Let B — UJ2 r U T . Then, the rank of the block symmetric matrix 



Z = 



A B 
B T A 



(2.8) 



equals the rank of A, i.e., rank(iT) = rank(A). 

Proof. First, recall two basic facts (see e.g.. lHorn and Johnson 1985| ) about ranks of matrix products: 

(i) Let X be m x n and Y any n x p matrix with rank n. Then, rank(Xy) = rank(X); and 

(ii) Let X be m x n and Y any p x m matrix with rank m. Then, rank(KX) = rank(X). 
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It is easy to see that Z may be factorized as 



Z 



UU T UJ r U T ' 
UJlU T UU T 







~u 


0" 


T>2r 




" 









u 








where Li r 



— J-lr 



Since U has rank 2r, the rank of the direct sum rank(?7 © la) = 4r; similarly rank({/ T © ^2) = 4r. Thus, using 
Properties (i) and (ii) of ranks, we conclude that rank(Z) = rank(L2 r ). But Since l2 r and Jir are invertible, 
elementary manipulations show that rank(P2r) = ranker) + ranker — J2rl2r(—J2r)) — ranker) = 2r, 
which follows upon noting J2 r J2r = —Izr- D 

As a consequence of [5] we immediately have the following corollary. 

Corollary 3. Let P be as defined by {2.$. Then, rank(P) = 2. 



Proof. Notice that fiy = cos(xi— Xj) = cos(xi) cos(a;j)+sin(xi) sin^), while 6^ = sin(xi—Xj) = sm(xi) cos(xj)- 
cos(xi) sm(xj). Thus, A is of the form UU T , where rank([/) = 2, while B is of the form UJ2U T , for a suitable 
nx2 real matrix U. Thus, rank(P) = rank(A) = 2. □ 

3 Explicit eigenvalues 

In this section we derive explicit expressions for eigenvalues of P. Here Corollary [3] plays a key role since it 
allows us to conclude that rank(P) < 4 (using (|2.7p ). Thereafter, basic linear algebra shows that the nonzero 
eigenvalues of P lie among the roots of its reduced order-4 characteristic polynomial: 



XP(X) := A 4 - 0i A 3 + 2 A 2 - 3 A + 4 , 



(3.1) 



where m is the mth elementary symmetric polynomial of P. Using (I3.1[) we can now proceed onto computing 
the eigenvalues of P. We will see that when rank(i) = 2, the characteristic polynomial (|3.1[) actually has 4 
nonzero roots, thus actually the equality rank(P) = 4 holds. 

All that remains is to factorize the polynomial (|3.1[) . which will then yield the desired eigenvalues. To that 
end, we first compute the coefficie nts m , a job that is simplified by invoking the well-known Newton's identities 



for symmetric functions [see e.g.. iMeacu . Il992l |. and the basic fact that ^Af i (P) = trP" 1 . Specifically, we 
have the following identities: 



01 = 


trP 






202 = 
303 = 


(trP) 2 - trP 2 

02 trP-trPtr P 2 


+ tr P 3 


(3.2) 


404 = 


03 trP-02 trP 2 + 


trPtrP 3 - 


trP 4 . 


(1 < m < 4) via ([3~2]) 


we require knowledj 


;e of tr(P m ) 


. And it is at this point where we may 



To compute 

hope for benign simplification due to the special trigonometric structure of P. Here are the details. 
Lemma 4 (Powers). Let P and L be as defined by equations \2.1)) and (|2.3[) . respectively. Then, 

tr(P m ) = 2 tr(L m ), m = 1, 2, 3, 4. (3.3) 

Proof. The case m — 1 is obvious, since an ~ 1. For m > 1, we propose computing the following formulae: 

tr(P 2 )=^. . p ijPij (3.4) 

tr ( p4 )=E^ =1 4-4- (3-6) 

We split each of traces above into 4 parts, corresponding to the 4 blocks of P. Pictorially, 
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p ^> 



Si 


s 2 


s 3 


s 4 



where Si (l < i < 4) denotes the contribution to the trace from the corresponding block of P" 
(1) The case m — 2. 



tr 



( P ) = >.. . . PijPij 
* — '»,.7 = 1 



n In 



2n n 



2n 2n 



E ^+ E E £4 



z,j — 1 i— 1 j—n-\-l 

= 5i + ^2 + £3 + S-4- 



i— n+1 j— 1 



£=71+1 j— n+1 



We compute each of the sums Si now. The algebraic simplifications exploit two main points: (i) symmetry of 
L; and (ii) the fact that sin(x) is an odd function. 

n n n 

& = £ 4 = £ a % = £ 

i)J=l * ) J=1 H?=l 

n 2n n n n n 

s 2 = £ £ 4=££4=££^ 

2—1 J =n+l 1=1 j — 1 1=1 J = l 

2n n n n n 

^3= E £4=££^=£^ 

i=n+l j — 1 i=l j=l ijj = l 

2n 2n n 

5 4 = £ = £ 4 = 5i - 

z=n+l j— n+1 i } j=X 

Notice, although bij = ~bji, their squares are the same. Thus, S2 — S3, whereby we obtain 

tr(P 2 ) - 2(5i + S 3 ) =2Y,..{a%+b%)ll = 2^..(cos 2 ( a;< - Xj ) + sm 2 ( Xi - Xj ))l* = 2tr(L 2 ). 
(2) The case m = 3. 

For clarity, we introduce the shorthand fty = Xi — Xj. 



' A 


D 




' A 


B 




B T 


A 




B T 


A 





A 2 + BB T AB + BA 
B T A + AB T B T B + A 2 



tr(P 3 ) = [P 2 ]i^« = 5! + 5 2 + 5 3 + 5 4 

Direct multiplication shows that P 2 is given by 

p 2 

Using (|3.7p . we now compute each of the partial sums Si through S4. 

n 

Si = ([^ 2 ]ij + [BB T ]ij)a,ijlij = ^~],„ a>jjkj(cHka>kj +bikbjk) 

= / J.., aijUjkkhj{Q>ika>kj + i>ikbjk) — / J . , , ^jj^fe ^fcj (cos ^fc cos flfcj — sinfe sin#jy) 



S2 = / [AB + />'. 1 — > hjkj(a,ikkkhjbkj + hkkkhjQ>kj) 

= E^ijfe bijkjkkhj(0'ikbkj + bikCLkj) = £._. fe hjkjkkhj (cos 0^ sin^- + sin 6*.^ cosOkj) 

= £^- fc hjkjkkhj Sm(6ik + Okj) = £^- fc ^ijkjkkhj — £,y ^ijkj[L 2 ]ij ■ 



(3.7) 



Observe that the derivation below depends on bij — —bji, b k i = —bik, and bkj = —bjk- 
S 3 = ^ , , [B T A + AB T }ijbjj = V? bjikj y2 (fc/wajy + a ik b ok ) 

IJ IJ K 

— / J, , bjilijlikhj(bkiQ>kj + ttikbjk) = / J., hjhjkkhj {hk&kj ~ ^ik^kj) 

= 7 ^ bjikjkklkj(sm9ik cosOkj + cos^/c sin^fej) = / bijkjhkhj sin(0^ + (9/^-) 

Similar manipulations show that S4 — S±. Thus, we have 

tr(P 3 ) = 2(5 X + S 3 ) = X,/"-, + ^i)«tf [i 2 ]« = 2 tr(L 3 ). (3-8) 

TTie case to = 4. 



Given the derivations above, one may safely guess that 

tr(P 4 ) = 2tr(L 4 ). 

Let us explicitly see why this guess is true. As before, we write 

tr(P 4 ) = y [P% [P% 3 =S 1 +S 2 + S 3 +S. 



4- 

Si = ^JA 2 ] y + [BB T ] %J f = (J2 k aika W + btkb i k ) 2 
= (y^jr- kkhj(a>ikO'kj — bikbkj) 

y^ fy^ ; _ hkhj {cos 9ik cosOkj - sm9 ik sm9 k] )j = y^, ,. Wfcj cos(6> ifc + 0jy 



Similarly, by now we routinely see that 



2i2 



and also that S3 = S2 and S4 = S%. Therefore, we finally have 

tr(P 4 ) = 2(S 1 + S 2 ) =2j2..$j + = 2 £..[£ 2 ]? j = 2tr(L 2 L 2 ) = 2tr(L 4 ). □ 

Using Lemma U we obtain (j) m as shown by Proposition [5] below. 

Proposition 5. //rank(L) = 2, then the nonzero eigenvalues of the generalized trigonometric matrix P given 
by (OOP are 

Ai,2 = 7 + 8, X 3 a =7-^ (3-9) 

and i/ie resi 0/ i/ie eigenvalues are zero. 
Proof. First, we introduce the quantities 

7m :=2tr(i m ), to = 1,2, 3, 4, (3.10) 
and then plug them into p.2[) to obtain 

4>i =7i; 02 = 5(7? -72); ^3 = gTi - 57271 + l73! 404 = <^)37i -0272 + 7173 ~74- (3.11) 

We must now compute roots of the quartic polynomial XpW given by (|3.1j) . To obtain its roots explicitly, let 
us further refine the values that 4> m can take, by computing 7 m = A™(L) explicitly. 
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Thus, we first compute the eigenvalues of L. Recall that L = hl T + lh T is at most a rank-2 matrix. Thus, 
its nonzero eigenvalues must lie among the roots of the reduced polynomial 

XlW :=A 2 -0i(L)A + 2 (L)=O, (3-12) 

where <f>i(L) (i — 1, 2) are elementary symmetric polynomials of matrix L. Now define the quantities 

j = l T h, 6=[(l T l)(h T h)] 1 / 2 . (3.13) 

Then, 0i (L) = tr L = 2 7 , while 2 (L) = ^((tr L) 2 - trL 2 ) = 7 2 - S 2 , since 

tri 2 = tr(lh T + hl T ){lh T + hl T ) = tr(lh T lh + lh T hl T + hl T lh T + hl T hl T ) = 2( 7 2 + S 2 ). 

Plugging in these values into (|3.12[) we immediately see that the quadratic factorizes as 

XL (X) - A 2 - 2 7 A + 7 2 - S 2 = ( 7 - 5 - A)( 7 + 5 - A). (3.14) 

Thus, the (nonzero, unless j = S) eigenvalues of L are 

X 1=J + S, A 2 = 7 -d (3.15) 

Using (|3.15j) and performing some algebra we obtain the following equations: 

71 = 2(Ai + A 2 ) = 4 7 (3.16a) 

72 = 2(A 2 + A 2 ) =4( 1 2 + 5 2 ) (3.16b) 

73 = 2(A? + A|) = 4 ( 7 3 + 3 7 <5 2 ) (3.16c) 

74 - 2(Af + \f) - 4 ( 7 4 + 6~f 2 5 2 + S 4 ) . (3.16d) 

Using p,16ap - ()3.16dp and the definitions (|3.1ip . we obtain <f) m as follows 

0i = 47, 2 = 6 7 2 - 2,5 2 (3.17a) 

03 = 4( 7 3 - 7 5 2 ), 04 - (7 2 -^ 2 ) 2 - (3.17b) 

Recall now the reduced characteristic polynomial for P 

Xp(X) = A 4 - X A 3 + 2 A 2 - 3 A + 04, (3.18) 

into which we substitute (|3.17ap - (|3.17bj) and simplify to obtain the factorization 

X p(\) = {\-J + Sf(\- 1 -Sf. (3.19) 

This immediately yields the desired roots 

Ai,2 = 7 + 8, A 3i4 = 7-5. □ 

Remark: Observe that 7 = 5 can hold if only if rank(L) = 1 (we ignore the trivial case of rank-0 L). In this 
case P has only two nonzero eigenvalues. 



As a corollary we immediately obtain a proof of Conjecture [TJ which was posed bv lZhang et al.l 2012] , 

Corollary 6. Let the matrix P{oj) be as defined by (jl.ip and (|1.2p . Then, it has 2 positive and 2 negative 
eigenvalues given by 



n 



. /4n 2 -6n + 2\ x n /4n 2 -6n + 2\ 
A + (P) = 3 » - 1 + V 3 . A_(P) = - [n - 1 - . (3.20) 
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Proof. From (11.21) we see that L — le T + el T , so that S = V nl T l = ^/n\\l\\. Moreover, Zj — (i — l)/2 for 
1 <i < n. Thus, the eigenvalues of P are given by the roots of (|3.19[) . which equal 

A + =7 + Vn\\l\\, A_ =7- Vn||*||, (3-21) 
both with multiplicity 2. Since Z, = (i — l)/2, we have 

7 = Fe = |n(n - 1), ||Z|| = ^± { 2n^-3n^ + n) = ^y^3|^2 

Plugging these values of 7 and ||Z|| into p.21[) and simplifying, we obtain p.20p . □ 

We conclude by mentioning that instead of trigonometric matrices, we can carry out a similar derivation 
for more general rank-2 matrices that satisfy the hypotheses of Lemma [2] But the details are laborious, so we 
leave them as an exercise for the interested reader. 
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